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Abstract 

This study examines the simultaneous effects of heat and mass transfer on the three-dimensional boundary layer flow of 
viscous fluid between two infinite parallel plates. Magnetohydrodynamic (MHD) and thermal radiation effects are present. 
The governing problems are first modeled and then solved by homotopy analysis method (HAM). Influence of several 
embedded parameters on the velocity, concentration and temperature fields are described. 
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Introduction 

The boundary layer flows in the presence of heat and mass 
transfer are important in the metallurgical processes involving the 
cooling of continuous strips. At present this topic has been studied 
extensively via different flow configurations and assumptions. Few 
representative investigations in this direction can be found in the 
refs. [1-5] In all these studies the research has been carried out for 
flows induced by moving surface in semi-infinite expanse of fluid. 
Very little is said about such flows in the bounded domains. 
Dinarvand and Rashidi [6] examined a reliable treatment of 
homotopy solution for two-dimensional viscous flow in a 
rectangular domain by considering two moving porous walls. 
Axisymmetric finite element solution of non-isothermal parallel- 
plate flow has been presented by Zhang and Olagunju [7]. 
Martinez et al. [8] discussed the viscous flow in a curved channel 
with flexible moving porous walls. Bifurcation of multimode flows 
of viscous fluid in a plane diverging channel is examined by 
Akulenko and Kumakshev [9]. Mahmood and Asif [10] studied 
generalized three-dimensional channel flow due to uniform 
stretching of plate. Recently, an incompressible moving boundary 
flows with the finite volume particle method has been examined by 
Nestor and Quinlan [11]. 

To our information, the MHD channel flow due to stretching 
surfaces with heat and mass transfer is not investigated yet in 
presence of thermal radiation. Needless to say that mass transfer is 
involved in the net movement of mass from one locality to another 
in the system and has abundant applications in chemical 
engineering. Specific applications include evaporation of water, 
the diffusion of chemical contamination, separation of chemicals in 
distillation procedure etc. For separation processes, thermody- 
namics determines the amount of separation, while mass transfer 
determines the rate at which the separation is possible. Further the 
heat transfer is a process that concerns with the exchange of 
thermal energy from one physical system to another. The 



radiation effects are prominent in any transparent medium for 
instance solid or fluid, but may also even occur across vacuum in 
the form of electromagnetic rays. The purpose of current 
investigation is to examine the simultaneous effects of heat and 
mass transfer on the magnetohydrodynamic (MHD) three- 
dimensional channel flow of viscous fluid in the presence of 
thermal radiation. Section 2 includes the mathematical formula- 
tion and analysis with graphical results whereas section 3 includes 
the concluding remarks. 

Mathematical Formulation 

Three-dimensional flow of an incompressible viscous fluid 
bounded by two infinite parallel plates (separated by distance h) is 
considered here. The motion in fluid is due to equal and opposite 
forces along the x- and y-axes so that the plate is stretched in both 
directions keeping the origin fixed. The lower plate is a highly 
elastic membrane situated at z = 0 and the upper plate is uniformly 
subjected to constant injection in the channel fixed at z = h. A 
constant magnetic field of strength Bq is applied. In addition the 
heat and mass transfer effects in the presence of thermal radiations 
are also present. The flow is governed by the following expressions 

du dv dw _ , 1N 

du ^ du ^ du l dp ^ ^d 2 u ^d 2 u ^d 2 u^ oB\ 
dx dy dz pdx dx 2 dy 2 dz 2 p 

dv ^ dv ^ dv l dp ^ ^d 2 v ^ d 2 v ^ d 2 v ^ oB\ 

dx dy dz pdy dx 2 dy 2 dz 2 p 
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R = 5.0. A = 0.5 = 0. M = 1.0 =Sc = Pr. Rd = 0.3 



R = 5.0, X = 0.5 = 0, M = 1.0 = Sc = Pr. Rd = 0.3 
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Figure 1. "h curve of / at the 15th order of approximation. 

doi:1 0.1 371 /journal. pone.0086695.g001 



Figure 3. "h curve of (j) at the 15th order of approximation. 

doi:10.1371/journal.pone.0086695.g003 



dw dw dw I dp .d 2 w d 2 w d 2 w. ... 



JC ^dC w dc _ D d c 

dx dy dz dz 2 



u dT +y dT +w dT _ a d 2 T ^\6a*T^d 2 T^ 
dx dy dz m dz 2 3k*pc p dz 2 



(5) 



(6) 



subject to the boundary conditions 
u = u w (x) = ax, v = v w (y) = by, w = 0, C=C W , T=T W , at z = 0, 



electrical conductivity, <7* the Stefan-Boltzman constant, k* the 
Rosseland mean absorption coefficient, C the concentration of 
species, D the coefficients of diffusing species, T the temperature 
and ot m the thermal diffusity. The fluid phase temperature 
differences within the flow are sufficiently small so that T 4 may 
be described as a linear function of temperature. Hence last term 
in Eq. (6) is obtained by expanding T 4 in Taylors series about the 
free-stream temperature Tq and neglecting higher order terms. 

The boundary conditions for the present situation can be 
written in such a way that C w denotes the concentration at the 
surface, Co is the concentration far away from the sheet, T w is the 
surface temperature and Tq is the temperature far away from the 
surface. 

Using [10,17] 

n= j^z, u = axf(rj), v = ayg(rj), w= - ah{f \rj) + g(rj)} , 



u = 0, v = 0, w=-V 0 , C=C n ,T=T 0 at z = h. 



(7) 



In above equations u, v and w are the velocities in the x, y and z 
directions, respectively, Vo the constant injection velocity at upper 
wall, p the pressure, v the kinematic viscosity, p the density, a the 



T w — To 



C w — Co 



(8) 



the continuity equation is automatically satisfied and Eqs. (2) — (7) 
give 
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Figure 2. ti curve of g at the 15th order of approximation. 
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Figure 4. "h curve of 9 at the 15th order of approximation. 

doi:10.1371/journal.pone.0086695.g004 
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Table 1 . Convergence of the HAM solutions for different order of approximations when R = 5, 1 = 0.25, p = 0.3, M = 1 , Pr = 1 = Sc 
and Rd = 0.3. 



order of approximation 


-/(0) 


-g"(0) 




~0'(0) 


1 


3.902083333 


0.6584325397 


1.104166667 


1.104166667 


5 


3.894097726 


0.6612752030 


1.101226749 


1.072402632 


10 


3.894057150 


0.6612623541 


1.101228897 


1.071850269 


15 


3.894057129 


0.6612623520 


1.101228896 


1.071851257 


20 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 


25 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 


30 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 


40 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 


50 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 


60 


3.894057129 


0.6612623521 


1.101228896 


1.071851290 
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number, Rd = k * k u is the thermal radiation and Pr = ah j cc m is 

the Prandtl number. Note that the constants a>0 and b>0. 

Now we solve Eqs. (8 — 11) along with the boundary conditions 
(12) by using homotopy analysis method. The details of homotopy 
analysis method are already included in our various published 
articles ([12-20]) and references there in) hence these are excluded 
here in order to save space. We choose auxiliary parameters ry, h^, 
h^ and h# for the functions f, g, <j) and 0 respectively. These 
auxiliary parameters play a vital role in adjusting the convergence 
of the obtained series solutions. We have plotted h-curves in Figs. 
(1—4) to obtain the permissible values of these auxiliary 
parameters. It is found that range for admissible values of hy, h^, 
h^ and h# are — 1.5 <hy, h^, h^, < — 0.3 and — 1.2<h^— 0.3 It is 
noticed that series solutions converge in the whole region of r] 
(0<rj<\) for hy, h^, h«^ and h e = 1.0. 

Table 1 signifies that how much order of approximations are 
required for a convergent solution. Obviously 15th order of 
approximations are sufficient for the analysis under consideration. 
In order to validate our results, we have given a comparative study 
of present results with the existing results. The results are in an 
excellent agreement (displayed in Table 2). Figs. (5-9) are 
displayed for the effects of different parameters on velocity, 
concentration and temperature fields. The effects of M on / and 
g have been depicted in the Figs. 5 and 6. It is seen from Fig. 5 
that in view of an increase in Hartman number M, the velocity 
f (ff) decreases in the vicinity of stretching sheet (0 <f] <0.4) 
whereas / (f]) increases away from the stretching sheet 



Table 2. Comparison of values of -f"(0) and -g"(0) for viscous fluid (M = 0) (refs.[10]) when Re = 5.0, /? = 0.5 and 1 = 0.5. 





Order of approximation 


Present results 




ref. [10] results 




-/(0) 


-g"(0) 


-/(0) 




5 


2.987025 


0.5037159 


2.98703 


0.503716 


10 


2.990508 


0.5037637 


2.99051 


0.503764 


15 


2.990524 


0.5037629 


2.99052 


0.503763 


20 


2.990525 


0.5037629 


2.99053 


0.503763 


25 


2.990525 


0.5037629 


2.99053 


0.503763 
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r-Reif'/ -if -ff -g/")-MY =0, (9) 
g iv -Re(g'g" -fg" -ff -gg"')-M 2 g"=0, (10) 

4>" + Sc[<f+g)4>']=0, (11) 
(l + *Rd)e" + Pr[(f+g)e']=0, (12) 

/ (0) = 1 , g (0) = ft, f + g = 0, 0(0) = 1 , 9(0) = 1 , 

/'(1) = 0, g'(l) = 0,f+g = A, 0(1) = O, 0(1) = 0, (13) 

where KQ = ah 2 /v is the Reynolds number, P = b/a(a^§) is the 
ratio of stretching coefficients, X=Vo/ah is the dimensionless 

injection parameter, M = u u , Sc = ah r/D is the Schmidt 
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A = 0.3 = p. Re = 5.0 




0 0.2 0.4 0.6 0.8 1 



Figure 5. Effect of M on the velocity component /. 

doi:1 0.1 371 /journal. pone.0086695.g005 

(0A<rj< 1). This is in accordance with the reason that magnetic 
field retards the fluid particles and slows down the motion in the 
vicinity of stretching sheets but obviously to satisfy mass 
conservation constraint a decrease in fluid velocity in the vicinity 
of stretching sheets is compensated by an increase in fluid velocity 
in the upper half of channel due to constant injection at the upper 
wall. This gives rise to a cross-over behavior which is obvious from 
Fig. 5. Variation of M on f and g are quite similar (see Fig. 6). 
The variation of Sc on the concentration field (j) is shown in Fig. 7. 
It is seen that concentration field (j) and concentration boundary 
layer are decreasing function of Sc. Since Schmidt number Sc is 
depending inversely on diffusion coefficient D. Hence an increase 
in Sc shows a decrease in diffusion D. This ultimately shows a 
decrease in concentration field <j). Fig. 8 elucidates the effects of 
Prandtl number Pr on the temperature field 6. Here temperature 
field 6 and thermal boundary layer thickness are decreasing 
function of Pr . This is because of the reason that larger Prandtl 
number corresponds to the weaker thermal diffusivity and thinner 
boundary layer. The influence of thermal radiation parameter Rd 
on 6 are given in Fig. 9. It is observed that thermal radiation 
causes a slight increase in temperature. 

Concluding Remarks 

The present study describes the effects of heat and mass transfer 
on the magnetohydrodynamic (MHD) three-dimensional bound- 



X = 0.3 = p, Re=5.0 
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Figure 6. Effect of M on the velocity component g . 
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X = 0.3 = /3, M= 1.0. Re=5.0 
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Figure 7. Effect of Sc on 0. 

doi:10.1371/journal.pone.0086695.g007 



X = 0.3 = p,M= 1.0. Re = 0.5. Rd=0.5 
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Figure 8. Effect of Pr on 6. 

doi:10.1371/journal.pone.0086695.g008 



X = 0.3=p,M= 1.0. Re = 5.0. Pr = 0.5 
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Figure 9. Effect of R d on 9. 

doi:10.1371/journal.pone.0086695.g009 
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ary layer flow of viscous fluid between two infinite parallel plates. 
The main points are summed up as follows: 

• It is observed that magnetic field M retards the flow near the 
stretching boundary. 

• Injection causes as increase in velocity to compensate the effects 
of magnetic field so that to mass conservation constraint is 
satisfied. 

• Prandtl number Pr and radiation parameter Rd show opposite 
behavior on temperature field 6. 

References 

1. Sahoo B, Do Y (2010) Effects of slip on sheet-driven flow and heat transfer of a 
third grade fluid past a stretching sheet. Int. Com. Heat Mass Transfer, 37: 
1064-1071. 

2. Jamil M, Fetecau G, Imran M (201 1) Unsteady helical flows of Oldroyd-B fluids. 
Comm. Nonlinear Sci. Numer. Simulat. 16: 1378-1386. 

3. Wang S, Tan WG (2011) Stability analysis of soret-driven double-diffusive 
convection of Maxwell fluid in a porous medium. Int. J. Heat Fluid Flow, 32: 
88-94. 

4. Hayat T, Ali S, Awais M, Obaidat S (2013) Stagnation point flow of Burgers' 
fluid over a stretching surface. Progress Gomp. Fluid Dyn. 13: 48—53. 

5. Awais M, Alsaedi A, Hayat T (2013) Unsteady flow of Maxwell fluid with 
thermal-diffusion and diffusion-thermo effects. Int. J. Numer. Meth. Heat Fluid 
Flow, Inpress. 

6. Dinarvand S, Rashidi MM (2010) A reliable treatment of a homotopy analysis 
method for two-dimensional viscous flow in a rectangular domain bounded by 
two moving porous walls. Nonlinear Analysis: Real World Appl. 1 1: 1502-1512. 

7. Zhang S, Olagunju DO (2005) Axisymmetric finite element solution of non- 
isothermal parallel-plate flow. Applied Math. Gomp., 171: 1081-1094. 

8. Martinez DM, Seo DJ, Jong J (2008) Viscous flow of a Newtonian fluid in a 
curved channel with flexible moving porous walls. Chemical Engng. Sci., 63: 
1882-1889. 

9. Akulenko L, Kumakshev S (2008) Bifurcation of multimode flows of a viscous 
fluid in a plane diverging channel. J. Appl. Math.Mech. 72: 296-302. 

10. Mehmood A, Ali A. (2007) Analytic homotopy solution of generalized three- 
dimensional channel flow due to uniform stretching of the plate. Acta Mech Sin, 
23: 503-510. 



• Schmidt number Sc causes a decrease in concentration field (j). 



Author Contributions 

Conceived and designed the experiments: MA AS. Analyzed the data: MA 
TH AS. Contributed reagents/materials/analysis tools: MA. Wrote the 
paper: MA. Provided financial support: AA. 



1 1. Nestor RM, Quinlan NJ (2010) Incompressible moving boundary flows with the 
finite volume particle method. Gomp. Methods Appl. Mech. Eng., 199: 2249 
2260. 

12. Liao SJ (2004) On the homotopy analysis method for nonlinear problems. Appl. 
Math. Gomp. 147: 499-513. 

13. Abbasbandy S, Shirzadi A (2011) A new application of the homotopy analysis 
method: Solving the Sturm Liouville problems. Comm. Nonlinear Sci. Num. 
Simulat. 16: 112-126. 

14. Abbasbandy S, Hayat T (2011) On series solution for unsteady boundary layer 
equations in a special third grade fluid, Comm. Nonlinear Sci. Numer. Simulat. 
16. 3140-3146. 

15. Bataineh AS, Noorani MSM, Hashim I (2008) Approximate analytical solutions 
of systems of PDEs by homotopy analysis method. Gomp. Math. Appl. 55: 
2913-2923. 

16. Bataineh AS, Noorani MSM, Hashim I (2009). On a new reliable modification 
of homotopy analysis method. Comm. Nonlinear Sci. Numer. Simulat. 14: 409 
423. 

17. Hayat T, Awais M, Obaidat S (2012) 3-D flow of ajeffery fluid in a channel of 
lower stretching wall. Euro. Phys. J. Plus, 127: 128-138. 

18. Hayat T, Awais M, Asghar S, Obaidat S (2012) Unsteady flow of third grade 
fluid with Soret and Dufour effects. ASME, J of Heat Transfer, 134: 062001. 

19. Rashidi MM, Pour SAM, Hayat T, Obaidat S (2012) Analytic approximate 
solutions for steady flow over a rotating disk in porous medium with heat transfer 
by homotopy analysis method. Computers & Fluids, 54: 1-9. 

20. Rashidi MM, Pour SAM, Abbasbandy S (201 1) Analytic approximate solutions 
for heat transfer of a micropolar fluid through a porous medium with radiation. 
Comm. Nonlinear Sci. Numer. Simulat., 16: 1874-1889. 



PLOS ONE | www.plosone.org 



5 



January 2014 | Volume 9 | Issue 1 | e86695 



